Improved exponential stability for 
near-integrable quasi-convex Hamiltonians 

Abed Bounemoura * and Jean-Pierre Marco ^ 
April 12, 2010 

Abstract 

In this article, we improve previous results on exponential stabil- 
ity for analytic and Gevrey perturbations of quasi-convex integrable 
Hamiltonian systems. In particular, this provides a sharper upper 
bound on the speed of Arnold diffusion which we believe to be opti- 
mal. 

1 Introduction 

This paper deals with some stability properties of near-integrable Hamilto- 
nian systems of the form 

(H{0,i) = hii) + fie,i) 

\\f\<e«l 

where {0,1) £ T" x are action- angle coordinates for the integrable part 
h and / is a small perturbation, of size e in some suitable topology defined 
by a norm | . | . If the system is analytic and under a suitable quantitative 
transversality condition called steepness, Nekhoroshev ( |Nek77j . |Nek79j ) 
proved that the action variables I{t) are stable for an exponentially long 
time, in the sense that for e sufficiently small, one has 

\I{t) - lol < cie\ \t\ < C2 expicae"") 

for any initial action Iq. The constants 01,02,03, a and b depend only h, 
one calls a and b the stability exponents and among them the value of a 
is of course the most important one, as it specifies the time-scale of sta- 
bility. Nekhoroshev's estimates complement the well-known KAM theory 
( |Kol54] ■ see also |Pos01| for a nice survey) which gives, under some mild 
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non-degeneracy condition on h and for e small enough, the existence of a 
constant c such that 

|/(t)-/o| <cVI 

for any time t G M, but only for a strict subset (of large relative measure) 
of the set of initial conditions. Of course KAM theory gives much more 
information, these stable solutions are in fact quasi-periodic and y^-close 
to the corresponding unperturbed solutions. In particular, for n = 2 and in 
the case when h is iso-energetically non-degenerate, this stability property 
even holds for all solutions. On the contrary, for n > 3, following Arnold 
( |Arn64| ) one can find examples of near-integrable Hamiltonian systems with 
a solution satisfying 

- /ol > 1 

for some time r = T(e) > 0, no matter how small the perturbation is. 
Such instability is commonly referred to as Arnold diffusion. Obviously 
Nekhoroshev's estimates give a lower bound on the diffusion time T (or 
equivalently an upper bound on the diffusion speed) which is exponentially 
large (or exponentially small when referring to the rate of diffusion). 

This paper is concerned with the precise time-scale at which stability 
breaks down and instability takes place, by which we mean the precise value 
of the exponent a, in the special case where the unperturbed Hamiltonian 
h is quasi-convex (that is convex when restricted to its level subsets). 

The quasi-convex case, which is of both practical and theoretical interest, 
has been widely studied in Nekhoroshev theory, essentially for two reasons. 
First the proof is much easier in this situation and a more refined result is 
available: the stability exponents may be chosen as 

a = b= (2n)~^ 

These facts are best illustrated by the striking proof given by Lochak ( [Loc92| . 
see also |LN92j and [LNN94J ) , though they are also accessible via a more tra- 
ditional approach as was shown by Poschel ( |Pos93] ). Note also that these 
values have been generalized by Niederman in the steep case ( |Nie04j ) using 
both ideas of Lochak and Poschel. Yet there is another reason for which the 
quasi-convex case is interesting, which is the so-called stabilization by reso- 
nances: if a solution starts close to a resonance of multiplicity m, m < n, 
then it posses better stability properties, described by the "local" exponents 

am = bm = (2(n - m))"\ 

This is a quite surprising fact, as it shows that even though resonances are 
the main cause of diffusion, at the same time they improve finite time stabil- 
ity. However, this property certainly does not hold without some convexity 
assumption (in the steep case for instance). 
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The optimality of the exponent a, in connection with the maximal speed 
of Arnold diffusion, has been first studied by Bessi who introduced power- 
ful variational methods to revisit Arnold's example and estimate the speed 
of diffusion. In |Bes96j and [Bes97j , he proved that the latter is of order 
exp (e~2) for n = 3 and exp (e for n = 4. Moreover, in Bessi's exam- 
ple the solution passes close to a double resonance, and so the speed is the 
highest possible in this case, in view of the values of the local exponent 02 
for n = 3 and n = 4. Recently, using similar variational arguments, this 
result has been generalized to an arbitrary number of degrees of freedom 
n by Ke Zhang ( |Zha09] ). namely he constructed a special orbit passing 
close to a double resonance for which the speed of diffusion is estimated by 
exp (e~2(^). 

Another approach has been proposed by Marco-Sauzin ( |MS02j ) and 
Lochak-Marco ( [LM05] ). following novel ideas of Herman. In |MS02] the 
authors showed that Nekhoroshev's estimates extend to perturbations of 
quasi-convex Hamiltonians which are a-Gevrey regular, a > 1, with the 
exponents 

a = (2an)~\ b = {2ny^ 

and local exponents 

am = {2a{n - m)y'^ , 6^ = (2(n - m))"\ 

Note that 1-Gevrey functions are exactly analytic functions, and basically 
when a ranges from one to infinity a-Gevrey functions interpolate between 
analytic and C°° functions. Therefore this result generalizes the estimates in 
the analytic case. Using a geometric mechanism different and more precise 
than Arnold's one, in |MS02| the authors constructed a drifting orbit with 
speed of order exp (e 2a(n-2) ^ jj^ the non-analytic case, that is when a > 1. 
Adding some more technical ideas, it was shown in |LM05j that the example 
also works in the analytic case but the speed was estimated as exp (e 2(n-3) ^ ^ 
which is only close to optimal (however refinements are certainly possible to 
reach the value (2(n — 2))~^ in this class of examples). 

Therefore, if the unperturbed Hamiltonian is quasi-convex, the best ex- 
ponent of stability a up to now satisfy 

{2ny^ <a< (2(n - 2))'^ 

in the analytic case and more generally 

(2an)~^ < a < {2a{n - 2))~^ 

in the Gevrey case. The goal of this paper is to improve the lower bound 
both in the analytic or Gevrey case, so as to have 

(2(n - 1))-^ -6<a< (2(n - 2))-^ 
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and 

(2a(n - 1))-^ - 5 <a< {2a{n - 2))"^ 

for S > but arbitrarily small (see Theorem 12.11 and Theorem 12.21 in the 
next section). We believe that this bound is optimal, in the sense that one 
could reach the value (2(n — 1)) in Arnold diffusion, using a significantly 
different mechanism of instability. 

2 Main results 

1. In order to state our main results, let us now describe our setting more 
precisely, beginning with the analytic case. Let B = B{0,R) be the open 
ball of M" of radius R > 0, with respect to the supremum norm, centered 
at the origin. Given s > 0, we let AsiV) the space of bounded real-analytic 
functions on P = T" x i3 which extend as holomorphic functions on the 
complex domain 

Vs = {{9,1) E (C^/Z") X I \I{9)\ < s, d{I,B) < s}, 

and which are continuous on the closure of Vg. Here we denoted by T{9) 
the imaginary part of 0, by | . | the supremum norm on C" and by d the 
associated distance on C". It is well-known that AsiT>) is a Banach space 
with its usual supremum norm | . |s, where 

I/I, = sup |/(z)|, feAs{V). 

In the following, we shall denote by 

5, = {/ e R" I d{I,B) < s} 

the real part of our domain Ds in action space. The geometric parameters 
n, R, s are assumed to be chosen once and for all in the following. 

We now introduce the parameters related to the choice of the system. 
The integrable part /i : i?, — )• M will be assumed to be strictly quasi-convex: 
the gradient map V/i does not vanish and there exists a positive number m 
such that 

V'^h{I)v.v > m\v\'^ (QC(m)) 

holds for any I £ Bg and any v orthogonal to V/i(I) (with respect to the 
Euclidean scalar product). Moreover, the derivatives up to order 3 of /i on 
Bs are assumed to be bounded: there exist M > such that for all I € Bg, 
one has 

|9^/i(/)| <M, 1 < \k\ < 3. {B{M)) 
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Therefore we will consider systems of the form 



H{9,I) = h{I) + f{9,I), 

(C{M,m,e)) 

/i satisfies (QC(m)) and (S(M)), 

J/|s<e. 

Note that we get rid of the geometric parameters in the notation. In the 
following we will call stable constant (in the analytic case) any positive con- 
stant c which depends on the whole set of parameters, that is n, R, s, M, m, 
together with a parameter 5 or /o to be defined below, but not on a particular 
choice of H satisfying condition (C(M, m,e)). 

2. The main result of the paper is the following. 

Theorem 2.1. Consider a real number 5 satisfying 

< 5 < (2n(n- 1))-^ 

Then there exist stable constants ci, C2, C3 and Eq such that if < e < £q, 
and if H satisfies {C{M,m,e)), the following estimates 

\I{t) - Io\ < cie^^'^~^\ \t\ < C2exp (^cse^^^^ 

hold true for every initial action Iq G B{0,R/2). 

Moreover, consider a real number p satisfying < p < R/2. Then there 
exist stable constants C4, C5 and eq such that ifO<e<eo, then 

- M < P, \t\ < C4 exp (c5e" 

for every Iq G B{0,R/2). 

Choosing our constant 5 arbitrarily close to zero, our result ensures sta- 
bility for a time-scale which is arbitrarily close to exp (e ^(n-i) ^ ^ therefore we 
improve the previous results of stability obtained independently by Lochak- 
Neishtadt f |Loc92j and [LN92]) and Poschel f |Pos93j ^. which were believed 
to be optimal. 

In fact in the extreme case where 6 = (2n(n — which in our 

situation gives the worst stability time (but of course the best radius of 
confinement), our result reads 



1 , , / J_N 

2n 



- M < cie2n , |t| < C2 exp (cse" 

and we recover the previous result of stability. Hence, when our parameter 
6 ranges from (2n(n — 1))~^ to zero, our Theorem "interpolates" between 
previous results of stability and what should be the optimal stability. 
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Indeed, in the other extreme case which corresponds to the second part 
of our Theorem, our result does not give stability since the radius of confine- 
ment can be arbitrarily small but no longer tends to with e. We believe 
that this is not an artefact of the method and that instability should oc- 
cur at this precise time-scale, at a time of order exp (e 2(n-i)^_ \Ye plan 
to construct an example with an unstable orbit which has a drift of order 
one during such an interval of time. This necessitates to use a more refined 
instability mechanism in the neighbourhood of double resonances, a topic 
which is also crucial in connection with the problem of genericity of Arnold 
diffusion. 

3. Our result also holds if the Hamiltonian is only Gevrey regular. Let us 
recall that given a > 1 and L > 0, a function H G C°°{T>) is (a, L)-Gevrey 
if, using the standard multi-index notation, we have 

\H\a,L= Yl L\'\''{nr'"\d'H\v<oo 

where \ -{v is the usual supremum norm for functions on D. The space of 
such functions, with the above norm, is a Banach space that we denote by 
Analytic functions are a particular case of Gevrey functions, as 
one can check that G^'^{V) =Al{V). 

Let us introduce the main condition on the Hamiltonian systems in the 
Gevrey case 

'Hi9j) = hii) + fie,i), 

h satisfies (QC(m)) and (-B(M)), 
J/U,L < £■ 

We now call stable constant (in the Gevrey case) any positive constant 
c which depends on the whole set of parameters, that is a, L, n, R, s, M, m, 
together with a parameter 6 ov p which will be defined below, but not on a 
particular choice of H satisfying condition {C{a, L, M,m,e)). 

4. Our second result is the following Gevrey version of Theorem 12.11 
Theorem 2.2. Consider a real number 5 satisfying 

0<6< {2an{n - 

Then there exist stable constants c'l, C2, Cg and e'q such that if < £ < e'q, 
and if H satisfies C{a,L,M,m,e), the following estimates 

\I{t) - Io\ < c[e^(^, \t\ < 4 exp (4e"2«(n-i)+'5'j ^ 



{C{a,L,M,m,e)) 
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hold true for every initial action Iq € B{0,R/2). 

Moreover, consider a real number p satisfying < p < R/2. Then there 
exist stable constants c'^, C5 and Eq such that if < e < e'q, then 

\I{t) - Io\ < P, \t\ < C4 eXp (C5e"2a(n-1)) 

for every Iq G B{0,R/2). 

The same remarks as above apply in the Gevrey case. In particular 5 
can be chosen arbitrarily close to zero and our result ensures stability for 
an interval of time which is arbitrarily close to exp (e 2a(n-i)^^ However, 
our radius of stability is worse than in the analytic case, so we do not fully 
recover the result obtained in | MS02] . but of course the time of stability is 
the most important issue. 

5. To avoid cumbersome expressions in the following, when there is no 
risk of confusion we will replace the stable constants with a dot. More 
precisely, an assertion of the form "there exists a stable constant c such 
that / < eg" will be simply replaced with "/ <• g" , when the context is 
clear. Such modifications will only concern assertions stating the existence 
of stable constants, and dealing with equalities or (strict or large) inequalities 
of well-defined functions. 

6. In the rest of the paper, as usual and without loss of generality, we will 
only consider solutions starting at time t = and evolving in positive time, 
and initial conditions will be denoted by (Iqi^o) = (-^(O)i ^(0))- 

7. Finally in the rest of this text all the norms will be denoted by | . | , this 
will always be the supremum norm for vectors and the induced norm for 
matrices, except for vectors in for which | . | will stand for the ^""^-norm. 



3 The analytic case 

In this section, the geometric constants n,s,R together with the constants 
m and M are fixed once and for all. 

The proof of Theorem 12.11 relies on two elementary facts. The first 
one, which we recalled in the introduction, is the stabilizing effect of reso- 
nances: on account of quasi-convexity, solutions of any Hamiltonian satisfy- 



ing (C(M, m, e) ) and starting sufficiently close to a resonance are stable for 
a longer interval of time. Of course this concerns only some special solutions, 
but we need to consider all of them. Our second remark is that to deal with 
the remaining ones, one can take advantage of the geometry of resonances in 
the integrable system to obtain a confinement result. More precisely, using 
the iso-energetic non-degeneracy implied by our quasi-convexity assumption, 
we will show that solutions that avoid all resonances are necessarily stable 
for all time. 



7 



1. Let us begin by making the first point explicit, using Poschel's approach 
of Nekhoroshev's theory. We shall denote by = V/i(i?) the space of 
frequencies. Let A be a sub-module of of rank r, with r £ {1, . . . ,n}. 
The resonant space associated with A is defined by 

Ra = {uj e n \ k.uj = 0, V/c G A}. 

Given a real number K > 1, we will say that A is a JC-sub- module if it 
admits a Z-basis {k^, . . . , k^} satisfying | A:*| < K for i G {1, . . . , r}, where 

1^ I ~ 1^1 1 ~l~ ■ ■ ■ ~l~ l^nl' 

As usual, it is enough to consider only maximal sub-modules, which are 
those that are not strictly contained in any other sub-module of the same 
rank. Given such a sub-module, we define its volume by 

|A| = Vdet *MM 

where M is any n x r matrix whose columns form a basis for A (this is 
easily seen to be independent of the choice of such a matrix) . The following 
stability Theorem is due to Poschel. 

Theorem 3.1 (Poschel). Let A be a K -sub-module of of rank r, with 
r G {0, . . . , n — 1}. Assume that e > and K > 1 satisfy 

and H satisfies ( fC(M, m, e)\ ). Then, if u;(0) = V/i(Io), for any solution 
{0{t), I{t)) with Iq G B{0, R/2) and d{u}{0), Ra) <• y/e the following estimates 

\Iit)-Io\< (e|A|2)5(^, i< exp(-e|A|2)"5(^^^, 

hold true. 

As we recalled in the introduction, we use a dot in the various inequalities 
to abbreviate an assertion such as "there exists a stable constant Cj such 
that", located at the beginning of the statement. 

The previous Theorem exactly gives the content of Theorem 3 in |Pos93j , 
to which we refer for a possible choice of stable constants. Note that Poschel 
uses a more quantitative version of quasi-convexity, namely that there exist 
two positive numbers / and m such that at least one of the inequalities 

\Vh{I).v\ > l\v\, V'^h{I)v.v > m\v\'^, 

holds for any I £ Bg and any v G M". It turns out that this notion is in fact 
equivalent to our condition (QC{m)). 
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We shall only need Poschel's result in the special case where the K- 
sub-module A has rank 1 and where the solution starts precisely on the 
associated resonant manifold. So we let 



Rk = {]Rk 



where the union is taken over all (maximal) -R'-sub-modules of rank 1, so it 
is the set of simply resonant frequencies of order K. Moreover, the "volume" 
of a rank-one sub- module |A| is nothing but the Euclidean norm of its (two) 
Z-generators, so one gets the trivial estimate 

1 < |Ap < K'^ 

from which we deduce the following lemma. 

Lemma 3.2. Let K he a K -sub-module ofTL^ of rank 1. Assume that e > 
and K > 1 satisfy 

ei^2" <• 1 (1) 



and H satisfies {C{M,m,e) ). Then for any solution (9{t),I{t)) such that 
Iq E S(0,i?/2) and tij(O) G Rk the following estimates 



\I{t) - Iq\ <• {eK^) , t <• exp {-eK^) ^(n-i) 

hold true. 



2. Let us now turn to our second remark, which is a simple geometric 
property of the integrable system based on the iso-energetic non-degeneracy. 
The latter condition is known to be implied by quasi-convexity, and it can 
be interpreted in various ways (see ( |Sev06j ) and references therein), but for 
subsequent arguments we will adopt the following form. 

Lemma 3.3. Suppose h satisfies ( fQC(m)[ ). Then the map 

■■ BxRt — > M X (M" \ {0}) 
(/,A) ^ {hiI),Xuj{I)). 

is a local diffeomorphism in the neighbourhood of any point (Iq, Aq) G i?xR+. 

In particular, there exist stable constants po and C such that if p < po 
then is a diffeomorphism in restriction to the ball B{{Iq, Aq), p) ^ BxWf , 
whose image contains the closed ball -B(^'/i(/o), C/?) C M x M". 

The second statement is a consequence of the first one (see |MS02j lemma 
3.17 for quantitative estimates on the stable constants po and C), and the 
first statement is well-known but we give a proof for convenience. 
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Proof. By the inverse function Theorem, it is enough to prove that d(^j^^ \g^^h 
is non-singular at any point (/o,Ao) € B x M+. Given u £ W, v £ W^, we 
easily compute 

d{ioM)'^h{v,u) = {uj{Io).v,uu;{Io) + AoV^/i(/o)v) 

and we need to show that this vector is non-zero if the vector {v, u) G W^^^ 
is non-zero. If either u}{Iq).v 7^ 0, in which case the first component is non- 
zero, or = and hence the second component is non-zero (since « 7^ and 
so ULj{Io) / 0), the statement is obvious. Otherwise, u}{Io).v = and f / 0, 
since h satisfies {QC{m)) this gives 

{uuj{Io) + XoV'^h{Io)v).v = XoV^h{Io)v.v 

> Xom\v\'^ 

and therefore uu]{Iq) + AoV^/i(/o)?^ 7^ 0. □ 

3. We can now make one step further in the dynamical consequences of 
the structure of simple resonances. We consider a Hamiltonian H satis- 
fying (C(M, m,e)). We will actually focus on those simple resonances for 
which the ratio of two frequencies becomes rational. The following elemen- 
tary lemma will allow us to deal with this simple case. 

Lemma 3.4. Let I he a closed interval of length I > contained in [—1, 1]. 
Then there exist a rational number p/q G /HQ satisfying 

\q\ + \p\ < {4V2)l-^. 

The exponent in comes from the use of Dirichlet's Theorem on the 
approximation of real numbers by rational ones, but this result is not nec- 
essary as in the sequel a trivial bound of order l^^ would be enough. 

Proof. Let us write I = [x — 1/2, x + 1/2] for some x S [—1, 1], and let q be 
the smallest integer larger than \/2/~2, that is 

V2r^ <q< V2l-^ + 1. 

By Dirichlet's Theorem there exists an integer p G Z such that 

\x — p/q\ < 

Since \/2l~2 < q, we have q^^ < 1/2 and so 

\x -p/q\ < 1/2 

which means that p/q G I. Moreover, as / C [— 1, 1], then \p\ < q and 

\q\ + \p\ < 2q. 
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Recall that q < \f2l 2 + 1, but as / C [—1,1] we have / < 2 so that 
1 < \f2l~'^ and therefore 

q < {2V2)l-h. 

This gives 

\q\ + \p\ < (4V2)r5 
which concludes the proof. □ 

The following result is our main lemma. It essentially says that a (long) 
drifting orbit has to cross a simple resonance, since all other orbits are stable 
on the interval of time over which they are defined. 

Lemma 3.5. Consider e > and K > 1 such that 

K-^<1, eK^<l. (2) 



Let H be a Hamiltonian satisfying (C{M,m, e) ), r € U {+00} and let 
{9{t),I{t)) be a solution defined on [0,r[ with Iq G B{0,R/2). If 

uj{t) ^ Rk, 0<t<T 

then the inequality 

\I{t)-Io\<K-^, 0<t<T 

holds true. 

Once again, the exponent in comes from lemma [33] and hence from 
Dirichlet's Theorem, but a bound of order K^^ would be sufficient for the 
final result. Let us also add the if r is the maximal time of existence of the 
solution within the initial domain T" x B, then our stability estimate easily 
ensures that r = +00 which in turn implies stability for all time (see the 
proof of Theorem 12. ip . 



Proof. We will make conditions ([2]) explicit and prove that 
\I{t) - Io\ < 32 C-^K-^, 0<t<T, 

when 

< {32)-^Cpo, eK^ < 16, 



where po and C are the stable constants of lemma 13.31 

We will argue by contradiction, so we assume that there exists a time t 
for which 

\I{i)-Io\ > 32 C'^K~^. 

Consider the curve 

a{t) = {I{t),\uj{t)\-^) £ BxR+ 



11 



and let /> = 32 C^^i^^^ Then 

t* = inf{t e [0,t[ I a{t) ^ 5(f7(0),p)} 

is well-defined as the above set contains t. Now p < pQ so we can apply 
lemma [3^31 : the restriction of to the open ball B{a{0),p) is a diffeo- 
morphism whose image contains the closed ball B{^h{cr{0)), 32 R-"^). Con- 
sidering a slightly larger ball over which remains a diffeomorphism, this 
easily implies that 

^h{a{t*))^B{^h{am,32 K-^) 

that is 

- (/i(Io), 1^(0)1-1^^(0)) I > 32 
Using the preservation of energy one has 

|/i(/(t*)) -/i(/o)| <2e<32i^-2 

so that necessarily 

\\u{e)\-^uj{t*) - |a;(0)|-ia;(0)| > 32 K'^. 
Therefore there exists an index z G {1, . . . , n} such that 



OJ^{t*) W,(0) 



> 32 K'"^ 



\u{t*)\ |a;(0)| 

and this estimate means that the image of the interval [0, t*] under the 
continuous function 

|a;(t)| ^ ^ '''J 

contains a non-trivial interval / of length / = 32 . Now we can apply 
lemma [23] to find a rational number p/g G Q in reduced form and a time 
t' e [0, t*] such that 

|a;(t')| " q ^'^^ 

with 

H + <4v/2(32K-2)-| =Js:. (4) 

But as |w(t')l = l'^j(*')l some j G {1, . . . and replacing p with — p if 
ojj{t') is negative, the equality ([3]) can be written as 

qui{t')-pujj{t')=Q. (5) 

Now let us write k' = qci — pej G Z, then from dH) and ([5]) we have 

k'.uj{t') = 0, |A;'| < K, 

and since the sub-module generated by k' is maximal as p and (7 are co- 
primes, we have found that uj{t') G Rk- This gives the desired contradiction. 

□ 
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4. We can finally pass to the proof of Theorem 12.11 
Proof of Theorem \2.1l We choose K of the form 

with suitable stable constants Kq and Eq so that conditions ^ and ([2]) are 
satisfied if 

0<e<eo, 0<7<(2n)-^ 

With this threshold and these bounds on 7 both lemma [3^2] and lemma [531 
can be applied. Let {9q,Iq) g T" x B{0,R/2) and T the maximal time of 
existence within x B{0,R) of the solution {9{t),I{t)) starting at (^Oj-^o)- 
We have to distinguish between two cases. 

In the first case, we assume that uj{t) € Nk for all t < T. Then we apply 
lemma E3] with r = T to get 

|/(t)-/o|<-e^^, 0<t<T. 

From this estimate we can deduce that the solution I{t) belongs to some 
compact ball around Iq which, taking eq small enough (this is possible since 
7 > 0), is included in B{0, R). Therefore this solution is defined for all time, 
that is T = +00, and so the previous estimate gives 

-/ol <-e^^, 0<t<+oo. 

In the second case, there exists a smallest time < t* < T such that 
u){t*) belongs to the set Rk- We apply once again lemma [331 with t = t* 
to have 

-/ol <-e^^, 0<t<t*. 

Again, taking £q small enough we can ensure that I{t*) G B{0, R/2), then 
we can apply lemma [32] to the solution It*{t) = I{t + t*), which initial 
frequency belongs to Rk, to obtain 

\It'{t) - /t*(0)| <• [eK^]^^ , < t <• exp {-eK^Y^^^ . 

Setting 

1-27 
~ 2(n - 1) 

this gives 

-lt.(0)| <:e''\ < t<- exp(-e-"^). 
Since t* > 0, we get in particular 

\I{t)-I{e)\ <-e"^ r < t<- exp(-e-"^). 
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and we conclude that 

\I{t) - Io\ <• (e^^ + e"^) <• max{e2T,e'^^} , < t<- exp(-e-"^). 
Now, using the condition < 7 < (2n)~^ one has 

(2n)"^ <a^ < (2(n - 1))"^ 

in particular 

7 < 

and hence 

Therefore for all solutions, the estimates 

\I{t)-Io\<e^, < t<- exp(-e-"^), 

hold true. 

Finally, to obtain our statement just set 

5 = 7(n-l)-i 

so that 

= (2(n-l))-i -5 

llGnCG 

\I{t) - lol <e'^^-^\ 0<t< exp 
all this provided that 

0<5 < {2n{n - 1))"^ 

This ends the proof of the first part of the statement. For the second part, 
choosing K in terms of p but independent of e, the proof is similar and even 
simpler, so we do not repeat the details. □ 

4 The Gevrey case 

In this section, we will prove Theorem 1 2. 2 1 In fact, it will be enough to have 
a version of lemma 13.21 in the Gevrey the geometric considerations 

of the previous section still apply with no changes. 

Here we shall use a result from |MS02j . which follows the method in- 
troduced of Lochak ( [Loc92] ) from which results of improved stability near 
resonances actually originates. In the latter approach, the notion of "order" 
of a resonance is more intrinsic, however it is also more difficult to compute. 

Let A be a sub-module of Z" of rank r, with r G {1, . . . , n}, and choose 
a basis {k^,...,k^'} S for A. We define the matrix L of size r xn 
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with integer entries whose rows are given by the vectors fc* = {k\, . . . , A;^), 
1 < i < r, that is 

/kl • • • k^ 



L 



e M,,„(Z). 



G Mr 



\k[ ■ ■ ■ fej^ 

Then it is an elementary result of linear algebra that there exists integers 
di, . . . , dj. G Z, satisfying the divisibility conditions di \ ... \dr, such that L is 
equivalent to the diagonal matrix 

M ••• 0^ 

A = 

\0 dr 

Therefore one can write 

L = BAA (6) 

for some matrices A G GL(n,7j) and B G GL{r,'L). 

The numbers di are called the invariant factors of the module, and for 
a maximal module one can show that these numbers are all equal to one. 
The above normal form result can be proved equivalently by elementary 
operations on rows and columns or using the structure of finitely generated 
modules over a principal domain. 

One can easily check that ^A sends the standard sub-module (which is 
the one generated by the first r vectors of the canonical basis of Z") to 
the sub-module A. So quantitative informations about the sub-module are 
encoded in those matrices A G GL{n^'L). 

Following Lochak, we define ca (resp. c^) as the minimal value of the 
norm l^""*^] (resp. of |^|) among all matrices A G GL{n,'L) satisfying the 
relation ([6]) (it is easy to see that those constants depend only on A and not 
on the choice of such a matrix). In the space M„(Z) we may choose the 
norm | . | induced by the usual supremum norm for vectors, which is nothing 
but the maximum of the sums of the absolute values of the elements in each 
row. 

With those definitions, one can state the following stability result in the 
Gevrey class. 

Theorem 4.1 (Marco-Sauzin). Let K he a K -sub-module ofU^ of rank r, 
with r G {0, . . . , n — 1}. Assume that e > satisfies 

5(n— r) , -I 3(n— r) /2(n—r) , 

and H satisfies |C(q;, L, M, m, e j\ ). Then for any solution {d{t),I{t)) with 
Iq G B(0, R/2) and d(u}{0), Ra) <• ^/e the following estimates 



\I{t) - /ol <• 4/^4e2(„-,.) ^ t <• exp (-c 

hold true. 



5(n — r) \ 2a{n-r) 
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This is exactly the addendum to Theorem A in |MS02] , to which we refer 
for a possible choice of stable constants. Note also that in |MS02j the names 
of these constants are a bit different: there stands for what we have 
called CA, and another constant called ca is introduced which is obviously 
equivalent to ours. 

It will be sufficient to restrict ourselves to the case where the sub-module 
is of rank one. Even in this case, the constants ca and are in general very 
difficult to compute, so we will use only the obvious estimate 

CA < |^"^|, 4 < 

for some suitable matrix A satisfying ^ for a given i^T-sub-module A. 

Proposition 4.2. Let A be a maximal K -sub-module of rank 1. Then we 
have the estimate 

CA < < K. 

Let us point out that these estimates are very rough, however it seems 
difficult to improve the exponent in K. 

Proof. Let k = {ki, . . . , kn) G Z" \ {0} be a generating vector of A, with 
I A; I < K. Since A is maximal, the components of k are relatively prime and 
the invariant factor of A is equal to one, so 

B = il), A = (10...0), 

and a matrix A £ GL{n,'L) satisfies L = B/\A as in ([6]) if and only if its 
first row is equal to k. 

1. To prove our estimates, it will be enough to show that one can choose A 
such that the £^-norm of each of its rows is bounded by K, that is |j4| < K. 
Indeed, assuming the existence of such a matrix A, one immediately gets 

A < \A\ < K. 

Moreover, since the determinant of A is ±1, A~^ is a matrix of cofactors, 
therefore the absolute value of each of its element is trivially bounded by 
(n — 1)\K^~^, which gives 

CA < \A-^\ < nlK^'-K 

So it remains to prove that one can construct such a matrix, and we will do 
this by induction on n > 1. 

2. Let us state first an elementary remark. Let x £ Z^\{0} and y £ Z^\{0}, 
with d = gcd{x,y). Then there exist u,v satisfying the equation 

ux + vy = d 
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with the estimates 

\v\ , , \x\ 

To see this, note that the existence of at least one solution uo,vo follows 
easily from euclidean division algorithm. Then obviously 

i,y i.^ 

u = UQ-k-, v = vo-k-, 
a a 

is also a solution for any G Z. Therefore choosing k properly we can find 
at least one solution n, v with 

I I ^ 1 
\u\< — — 1. 

a 



For this solution one has 



since d—\x\< 0, so 



< d+ \ux\ 

< d+(M-l)|,:| 

< ¥w 



which proves our claim. Finally, if for instance y = 0, then d = x and one 
can obviously choose u = 1 and v = 0. 

3. Let i^T > 1 be given. We can now state our induction hypothesis H(n) 
for n > 1. 

H(n). Let k = {ki, . . . , kn) he vector of TI^ \ {0} with co-prime compo- 
nents such that \k\ < K. Then there exists a matrix A G GL{n, Z) with first 
row equal to k, which satisfies \A\ < K. 

The assertion H(l) is immediate since in this case k = (±1). Now for 
n > 2, assume that H(n-l) holds true and consider k = {ki, . . . ,kn) in 
Z"' \ {0} with co-prime components and |A;| < K. 

We may suppose that fc* = (ki, . . . ,kn-i) is non-zero (otherwise we 
consider fc* = {k2, ■ ■ ■ , kn)) and we set d = gcd(A;i, . . . , So c? > 1, the 

integers d~^ki,. . . , c/~^A;„_i are co-prime and 

\k.\ < -J. 
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By H(n-l) we can find a matrix 



hi 
\^n-l,l 



sucli tliat 



2,n-l 



^ra— 1,71— 1 / 



n-1 



G GL(n- 1,Z), 



i=i 

for each z G {2, . . . , n — 1}. 

Now since d and A;„ are co-prime, one can find integers u,v & 1^ such 
that 

ud + f A:„ = 1 

and therefore define a matrix 

^ ^1 ■ ■ ■ fen— 1 fen 

^2,1 • • • ^2,n-l 



A{u,v) 



^n-1,1 



^n— l,n— 1 
\n-l„,j-lu . I 1^n-l„ 



\(-l)"^S(i-^fel ••• {-lY'^vd-^kn-l {-lY-^uJ 

Expanding the determinant relatively to the last column easily proves that 
A{u,v) G GL(n,Z). 

As for the estimates, first assume that fe„ = 0. Then d = 1 and we may 
choose u = 1 and = 0, so obviously 

|A(1,0)| <i^. 

If now kn is non-zero, then by our previous remark we can choose 
so that 

^ I fen 1 1 1^*1 — 

which proves that the ^^-norm of the last row is bounded by K, and therefore 
[^(ti*,!;*)! < K. This ends the proof. □ 

With these estimates, one can deduce from Theorem 14.11 the following 
lemma. 

Lemma 4.3. Let A be a K -sub-module ofZ" of rank 1. Assume that e > 
and K > 1 satisfy 



and H satisfies (C{a,L,M,m,e)). Then for any solution {0{t),I{t)) such 
that Iq G B{0,R/2) and a;(0) G Rk the following estimates 

_/o|<.(eK("-i)(3n-i))2(^^, t<. exp(-ei^5("-i)'V^"'"-^) 

hold true. 
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The proof of Theorem 12.21 is now completely similar to the proof of 
Theorem 12. H but we shall repeat some details. 

Proof of Theorem \2.'A We choose K of the form 

with suitable stable constants Kq and Eq so that conditions ([2]) and ([7]) are 
satisfied if 

0<e<eo, < 7 < 5"^(n- 1)"^ 

Then we can apply both lemma 13.21 and lemma 13.51 in the same way as 
in the proof of Theorem 12.21 and setting 

_ 1 - 57(n - 1)^ _ 1 -7(n- l)(3n- 1) 
~ 2a{n - 1) ' ^ " 2(n - 1) ' 

we find that all solutions {6{t),I(t)) starting at {6o,Io), with Iq G B{0, R/2), 
satisfy 

\I{t) - Jol <• maxje^^^gfeTj ^ t< expi-e'''-'). 
Next using our condition 

< 7 < 5"^(n- 1)^2 

we have the bounds 

5~^(n - 2)(n - 1)-^ < < (2(n - 1))"^ 
hence 7 < 6^ so that 

\I{t) - Io\<e^', t< exp{-£-''-<). 

To conclude, just set 
so that 

= (2a(n - 1))"^ - 6 

tlGIlCG 

\I{t) - Io\ <• e-^'^^-'^-\ \t\ <• exp (^.e-(Mn-l))-'+5^^ 
all this provided that 

0<6< (2(n- 1))-^ 
from which we will only retain 

0<6< {2an{n - 1))"^ 

This concludes the proof of the first part, and the arguments for the second 
part are analogous choosing K that depends on p but independent of e. □ 
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